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Lorentz-Dirac equation is widely accepted as the classical equation of motion of 
an elementary point charge interacting with its own radiation (see for instance 



where = -F^v v is the external electromagnetic force. 

It is also well known that this equation is affected by some irreconciliable difficulties, 
that already show up in the case of rectilinear motion. Consider a free point charge that 
enters perpendicularly a parallel-plate capacitor at r = (proper time) and leaves it 
at 7~i > 0. For r < the charge is free, / M = and the solution to ^ is a uniform 
rectilinear motion, a M = 0. We can therefore take a M (0) = and 1^(0) = as 
initial data to integrate equation JH), so obtaining a unique solution for the velocity ?A 
Nevertheless, this solution has the drawback that, not only a M (r) does not vanish for 
r > n (when the external action has ceased), but it grows exponentially for r — > oo, 
what is known as runaway solution. 

Rohrlich [TH] put forward a way out consisting in that JTJ) is not the equation of 
motion, but it must be supplemented with an asymptotic condition: if the external 
force / M asymptotically vanishes, then the acceleration asymptotically vanishes too. 
As a result the resulting equation of motion is of integro-differential type and runaway 
solutions are ruled out (see also [T2]). 

This alternative however implies what is called preacceleration. Although the 
external force vanishes for r < 0, the solution to the above integro-differential equation 
presents non-vanishing acceleration before the force starts. This is not a surprising 
feature because, as pointed out in 0, it is a consequence of demanding the asymptotic 
condition in the future: the integro-differential equation of motion itself "foresees" what 
will happen in the future, r > T\. 

It thus seems as though we were facing the following dilemma pj: either (a) classical 
electrodynamics is self-contradictory or (b) Lorentz-Dirac equation is not the right 
equation that follows from classical electrodynamics. 

In view of this dilemma different stances are found in the literature. Rohrlich |IB| 
adopts the alternative (a) and adds that this is not a major trouble because the time 



dCHETlEni): 
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scale at which preacceleration shows up is too small (r ~ 10~ 23 s for electrons) far 
beyond the limits of validity of the classical theory. He further stresses that [20] «the 
notion of "classical point charge" is an oxymoron ...» since classical physics ceases to be 
valid below Compton wavelength. Moniz and Sharp also argued jTHlIIElIIZj that classical 
electrodynamics is only consistent in describing the motion of charges with radius larger 
than the classical electron radius, while the quantum theory of nonrelativistic charges 
is free of runaways and preacceleration. 

Other authors [23 Qlj embrace the alternative (b) on the basis that the derivation 
of Lorentz-Dirac equation involves Taylor expansions and therefore presumes that both 
the charge worldline and the external force are analytic functions. As a consequence, 
equation (JTJ is not valid in those points where x M (r) and / M (r) are not analytic. 
Particularly, Yaghjian [H3j studies a charged spherical shell of radius e and obtains 
an alternative equation: 



where T](r) = for r < and r/(r) = 1 for r > 2e/c. In another approach 
[EH EU El ED EH El EI|, the Lorentz-Dirac equation is thought of as a necessary — but 
not sufficient — condition the true equation of motion must fulfill. The true equation of 
motion, which will not have neither preacceleration nor runaway solutions, is of second 
order and can only be constructed by using a series expansion or a method of successive 
approximations. 

Others [2H] consider that the commented difficulties with Lorentz Dirac equation 
are not real physical problems, as they accept that acceleration can have a singularity 
in points where the applied force has a discontinuity. 

None of these justifications is fully satisfactory to us. Consider a classical charge 
modelled by a charge distribution and the corresponding energy-momentum distribution 
inside a sphere of radius e. Provided that a suitable set of constitutive relations is 
added, the local conservation of energy-momentum yields an evolution law for this 
continuous medium, which is deterministic and causal: the electric current and the 
energy-momentum distribution at t = determine the future values of these magnitudes. 
It is, to say the least, startling that, on taking the limit e — > 0, the causal and 
deterministic nature of the classical problem is lost. 
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Apparently Lorentz-Dirac equation is an unavoidable and flawless consequence of 
classical electrodynamics plus the local conservation of total energy and momentum 
[3HE1EB1- However, as the electromagnetic field contribution to the energy-momentum 
tensor is singular on the charge's worldline — it behaves as 0^ ~ 0(r~ 4 ) — some creative 
"tricks" are necessary to appropriately handle such a singular behavior in the energy- 
momentum balance. In our opinion, in most approaches to this problem some additional 
assumption slips into the reasoning through one of these "tricks". 

In this context, it is worth mentioning Rowe's work |211 [22], where more 
elaborated mathematical tools, namely regularization of generalized functions, are used 
to properly handle the singularity in G MI/ and obtain the Lorentz-Dirac equation. The 
use of generalized functions (or distributions) has also the advantage that no mass 
renormalization is necessary. 

We shall here use these same mathematical tools to review the derivation of 
Lorentz-Dirac equation and see that, contrary to the common belief, it is not a straight 
consequence of classical electrodynamics plus energy-momentum conservation, but it 
includes an elementary extra assumption. 

We shall here describe a point charge as a current distribution in an extended 
material body in the limit where the radius e — > 0. The total energy-momentum tensor 
results from two contributions: the electromagnetic part, 6^, which is associated to the 
field and pervades spacetime, and the material part, K^, which we assume confined to 
a world-tube of radius e and accounts for kinetic energy and the stresses that balance 
the electric repulsion among the parts of a neat total charge confined in a small volume. 

For e > both contributions 0^ and K^ v are continuous functions and can be 
considered separately. But in the limit e — > 0, the electromagnetic part presents a 
singularity 0(r -4 ) on the worldline. Therefore, in the limit e — > none of these two 
contributions can be properly defined, even resorting to generalized functions. However, 
nothing forbids the total energy-momentum tensor to converge to a generalized function 
for e — > 0, which will likely include 5 functions and its derivatives on the point charge 
worldline. 

In our approach we do not need to assume that the involved functions are analytic. 
Although Taylor expansions to some finite order are used, these hold for functions that 
are smooth enough, without need of analyticity [1]. 
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We shall examine what restrictions on the charge's motion follow from local 
conservation of energy and momentum, and find that the result is not Lorentz-Dirac 
equation but a somewhat less restrictive equation, already derived by Honig and Szamosi 
[TT] by extending Dirac's work. Then we shall see that this equation admits solutions 
that are free of both preacceleration or runaways. 

2. Statement of the problem 

2.1. Notation 

The retarded Lienard-Wiechert field of a point charge has an outstanding role along 
the present paper. Therefore it will be helpful to use retarded optical coordinates |2E] 
(as in ref. (21]) based on a timelike worldline Y = {^(r)} and an orthonormal tetrad 
{e^}o, = i 2,3,4, which is Fermi- Walker transported along T, 

= [ v »a v - v v a»\ e\ a) . (2) 

With a properly chosen initial tetrad, the latter evolution equation is consistent with 
the conditions 

e ( a fh^ = , <4) =v^ = ^ and a" = t>" , (3) 
where a 'dot' means «derivative with respect to r» and r]^ = (+ + H — ). Moreover, 
from now on we use units such that c = 1. 

For any point x in spacetime, the equation 

[x^-^(r)][^-^(r)]^ = 0, (4) 

supplemented with x 4 > z 4 (t), has always a unique solution, r = r(x), which defines a 
time coordinate for x. 

The space coordinates are 

X* = ef j} (x M - ^[r(x)]) (5) 

and the inverse coordinate transformation then reads 



^ - z^(r) + pv^r) + X^ {t) (r) , (6) 



where p = \\X\\ = ^(X 1 ) 2 + (X 2 ) 2 + {X 



3\2 
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The following relations and quantities, introduced in ref. [2E|, will be useful hereon: 
P = -K-^(r)K(r), F:=~[^-^(t)], 1 
n M := A; M — t> M , n M n M = 1 , fc^u^ = — 1 , J 

d^p = + p(a a n a )k^. (8) 

The unit space vector can be written as 

vi 

" _ p e (0 - U e (i) ■ 

Finally, the relationship between the volume elements in Lorentzian and in retarded 
optical coordinates is 

d A x = dr d 3 X = p 2 dr dp d 2 VL(h) , (9) 

where d 2 Vt{n) is the solid angle element. 

2.2. Some definitions and postulates 

A point charge is described by a current density four- vector, j^, and an energy- 
momentum tensor, , fulfilling 

d^f = , d^ v = and t"" = t u » , (10) 

respectively, the local conservation laws for total electric charge, energy-momentum, 
and angular momentum. 

We expect to obtain j M and t^ v as the limit of continuous distributions of charge 
and energy-momentum when the radius goes to zero, namely, 

(a) an electric current vector J M (e;x), which is confined to an "optical tube" of radius 
e around a timelike worldline T, that is, 

p(x) > e =>■ J^(e;x) = 0, (11) 

where p(x) is given by (JZJ), 

(b) an energy-momentum tensor T Ml/ (e;x) which results from two contributions: 

T^{e-x) = Q^ u {e-x) + K^{e-x). (12) 
The first term comes from the total electromagnetic field: 

F^(e;x) = F^(e;x)+FZ(x), (13) 
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namely, the sum of the retarded solution of the Maxwell equations for the current 
J M (e;x) plus an external free electromagnetic field. The second term in (fT2|) comes 
from the matter distribution which is also confined to the above mentioned "optical 
tube": 

p(x) > K^ u (e;x) =0. (14) 

The above continuous distributions of electric current and energy-momentum are 
assumed to fulfill the local conservation laws: 

8^ = 0, «9 M T^ = 0, T tu, = T v ' t . (15) 

We shall assume that both J M (e;x) and K fiv (e;x) are locally summable in R 4 and that 
F~dx) is continuous in R 4 . 



The retarded electromagnetic field is given by [HU 

F^(e;x) = ^ J J [l/ (e;x)d^D R (x-x')d A x' (16) 



with 



D R (x) = ^Y(x i )5(x p x p ) 



2tt" v " p> 

[y(x 4 ) is the Heaviside step function.] The retarded electromagnetic field is thus a 
continuous function and therefore locally summable in R 4 . 

In its turn, the electromagnetic contribution to the energy-momentum tensor, 



e^e;x) = l 



F^(e; x) (e; x) - -^F" a {e- x) F pa (e; x) 



(17) 



is also locally summable. 

The framework where the limits for e — > of J^(e,x) and T M ^(e;x) are 
mathematically meaningful and can be appropriately handled is the space £>'(R 4 ) of 
generalized functions (213 El- As locally summable functions, J M (e;x) and T M1/ (e;x) can 
be associated to generalized functions and, provided that the limits 

f = lim J"(e) G V'(R 4 ) , iT = limT^(e) G D'(R 4 ) 

exist, the continuity of differentiation operators in £>'(R 4 ) [HOI guarantees the 
conservation laws (TlT)|) as the limit of (fTH| for e — > 0. 

These conservation laws must now be understood in the sense of £>'(R 4 ), i. e. 
Vy? G P(R 4 ) , 

(d,f,<p) = and (8^,^ = 
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or 



and 



(lT ) d li <p) = Q. 



(18) 



3. The point charge limit 

3.1. The electric current 

If the support of J M (e;x) is the "optical tube" p(x) < e, then for any cp E £>(1R 4 ) such 
that suppy) does not intersect the worldline T, it exists E\ > such that ip(x) = 
whenever p(x) < e±. Therefore, for all e < ei, 



and in the limit e — > it follows that 

(J'*,(p) = 0, Vy?G£>(M 4 ) such that Tnsupp^ = 0. 

The support of the generalized function j M is therefore confined to the worldline T 
and, according to a well known result on generalized functions [31], can be written 
as a sum of (^-functions and its derivatives up to a finite order: 



with Z(«i-«»0f Vai = ; r = 1, . . . n. 

To model a point charge we only keep the lowest order term and, as a consequence 
of the conservation law (fTO . we have [26] 



where e is the electric charge of the particle and is a constant scalar. 
3.2. The energy-momentum tensor 

In our approach, the limits for K^ u (e) and B Ml/ (e) do not need to exist separately in 
T?'(1R 4 ). Our assumption is weaker and only the joint limit is assumed to be physically 
meaningful: 





J dr [P(t) S(x - z(t)) + P M (r) d a 5(x - z{r)) + ... 
+l ai - a ^(r)d ai ... a J(x-z(r))) 



(19) 




(20) 



1T = lim [K^{e) + 6^(e)] G £>'(K 4 ) . 



(21) 
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This fact expresses the notion that, although in the separate limits for both K fJ,v (e) and 
/iI/ (e) some infinities on the worldline Y could arise, these infinities will cancel each 
other, so that is defined in V'(R 4 ). 

3.2.1. The matter contribution If we restrict to test functions <p G V(W 4 — T), we have 
that 

limine) =0 GD'(l 4 -r). (22) 

Indeed, for any ip G V{R 4 — T) it exists €\ > such that (p(x) = whenever p{x) < e\. 
The confinement condition (fTljl then implies that 



Ve<e X) (K" v (e),(p)= / (Tx K^{e\ x) <p(x) = 
and equation (|22|) follows [32j. 



3.2.2. The electromagnetic contribution Recall now equations (fT6|) and (fTTj) . We have 
the pointwise limit 

limr( e; x) = F;(,) + F::(,), (23) 

where F^(x) is the retarded Lienard-Wiechert field, and is defined whenever x ^ T. It 
can be written as the sum of the radiation field plus the velocity field: 

Fj?{x) = Fr(x)+F%{x), (24) 

where, in the notation introduced in subsection 12.11 (also in ref. [26J): 

Ff{x) = j [(ah) vH^ + aH^] , (25) 

F%(x) = ^vW. (26) 

(Here (ak) = a x k\.) Similarly, for the electromagnetic energy-momentum tensor we 
have the pointwise convergence: 

lim e""(e;x) = Q^tx) , 

except at the points x G T. 

As a consequence of (J2BJ), Q^ u (x) can be splitted as 

e^(x) = e^(x) + bz( x ) + e^ x (x) . (27) 
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The first and second terms in the r. h. s. respectively result from substituting F R {x) 
and F^(x) into the quadratic expression ifTTjl. whereas 8^(1) comes from the cross 
terms. 

®rni x ( x ) an( i @ext( x ) are locally summable in M 4 . This is obvious for because 
it is continuous everywhere. As for Q^ ix (x), it is a sum of products of Fj££(x), which is 
continuous, and F R llJ (x), which is also continuous except for a singularity of order p~ 2 
on T that is cancelled by the factor p 2 in the volume element 0. Therefore, Q^ ix (x) is 
also locally summable in M 4 . We shall respectively denote: 

Ct := lim Oext(e; and 8^ ix := lim 6^ x (e; ar) (28) 

e— »0 e— >0 

with e ext , e mix e v(R 4 ). 

Let us now consider the Q R (x) contribution. It can be written as E 



6 «» = 

„2 



+ 

a^r + a u k^ - {an) {n^k u + n u k^)] + 



Airp 3 

4^ [a 2 - (an) 2 ] ^ , (29) 

which is continuous for x £ T. 

Owing to the p~ A and p~ 3 singularities on the r.h.s. of the above expression, not 
only Q R u (x) has a singularity on T, but in addition it is not locally summable. Therefore, 
no generalized function in V'{R 4 ) can be associated to Q R (x) in the standard way. 

Now, since Q R v (x) is a continuous function on M 4 — T, it is locally summable there, 
and this allows to take its finite part 9% G £>'(M 4 ) [33119]: 

(9%, <p) = J d 4 x &£{x) [<p{x) - Y(L - p) [<p{z) + pkfdMz)]] (30) 

for any (p G P'(M 4 ), where L is an arbitrary chosen length scale, z = z(t(x)) and r(x), 
k a and p(x) are defined in (jZj). 

Some points concerning the definition (|30j) are worth to comment: 

(i) The integral in the r.h.s. converges. Indeed, on the one hand, for p > L, Q R u (x) is 
continuous and ip(x) has compact support and, on the other, inside p < L we can 
use the mean value Taylor theorem [3] for the smooth function (p: 

(p(x) = <p(z) + pk x d x <p(z) + -p 2 k x k»d Xfl <p(z + p'k) 
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with < p' < p(x). Now, since •p is smooth and has compact support, dx^ip is 
bounded and it exists M > such that 

\ip(x) - [<p(z) + pk x d\(p(z)}\ < Mp 2 , xGsuppv?, < p < L . 

Hence the integrand in the r.h.s. of ijHUj) presents a singularity of order p~ 2 on V 
and therefore the integral converges. 

(ii) For a test function cp G X?(1R 4 — T), the function and all its derivatives vanish on T. 
Hence, (|30|) amounts to 



(iii) The definition (|3*Uj) consists of eliminating from the integrand as many terms in 
the Taylor expansion of <p(x) as necessary, in such a way that the remainder is 
summable and the condition (ii) above is fulfilled. As a consequence, the finite part 
9r G D'(R 4 ) is not unique. Indeed, on the one hand, we could have substracted 
some more terms in the Taylor expansion of (p, and obtained a convergent integral 
also fulfilling the requierement (ii). Besides, the length scale L is quite arbitrary 
and could even depend on r(x). 

This results in that 6^ is determined up to a finite sum of T-supported 8- 
functions and their derivatives, multiplied by arbitrary r-dependent coefficients, 
in an expression similar to (fTTIj) . We shall see that this lack of uniqueness in the 
definition of 9^ is not relevant at all, because we are not actually interested in 
9ft but in the total energy-momentum t^ . Here lies the difference between our 
approach and that of Rowe fl^\ . 

To give a more specific expression for 9 R V , we realise that since the r.h.s. of (|30|) is 
convergent, we can write 




(31) 



{0%, cp) = lim / d 4 xY(p - e) G%(x)<p(x) 




d 4 xY(p - e) Y(L - p) &%(x) [<p(z) + pk a d a ^{z)} 




which after a short calculation leads to 



0% = 0%- dr ([V^ - U^) 6{x - z(t)) - U x » v d x 5(x - z(r)j) , (32) 



where 




5(x - z) 



(33) 
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and the coefficients V^ v , U^ u and U Xflu depend on r and are: 

= - $- (3v"v u + rT) + ^ (5a W + 2a 2 rf v - a"a u ) , (34) 
6L 15 

= \e 2 L (a'V + aW) + (5a 2 */V + 2a 2 rT ~ ^a u ) , (35) 

U x » v = (3a V + 3a^ v - 2a V) + 

15 

2 r 2 

^— (2a 2 V + v v t^\ - a A [a'V + a V]) . (36) 

1 

Notice that they depend on the length scale L. 
We shall hereafter write 

= 9% + Ct + Cx • (37) 

Notice that 9^ G V'(R 4 ) C £>'(K 4 - T) . Now, since 0^(x) is locally summable in 
M. 4 — T, it can be considered as a generalized function 0^" G P'(M 4 — T) and, as a 
consequence of (pTT|) we have that 

0"" = 0^ in £>'(IR 4 - T) . 



3.2.3. The total energy-momentum tensor The total energy-momentum tensor is 
defined by the limit pTjl . For any test function (p G D(IR 4 — T) we have, as a consequence 
of (22) , that 

(t^, = lim y d 4 x 0^(e, x) 
and, using (j2*TJl . (f3Tf and (pTjl. we obtain 

Qr = v^p(t 4 -r). 

Therefore, t MZ/ — 6"^ G D'(IR 4 ) has support on T and, according to a well known result 
iU, it can be written as a finite sum: 



iT -9^ = j dr [wT{t) 5(x - z{t)) + m a ^(r) d a 5(x - z{r)) +... 

+m ai.. V(r) an< j (x _ z(r)) ] ; (38) 

where 

m {ai - ar)flu v ai =0, r = l...n. 

So far there is no correspondence between t^ u — 9^ v and the, so to speak, "matter 
contribution" to the energy and momentum. Therefore, we are not obliged to assign 
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this difference the value m / dr v^v v ' 5(x — z(t)), as it is done in ref. [2UE2T However, 
for the sake of the "elementarity" of the point charge we shall retain as few terms in (|38|) 
as possible, namely, 

*» = <r + / iT M - , W ) + m - W - , W) ] , 

which combined with l|32J) and (|3Tj) leads to 

^ = §<*> + Ct + Cx + tf f (39) 

with 

tf = y dr [tT{t) 5(x - z{r)) + p A ^(r) <9 A 5(x - z(r))] (40) 

and 

p^ = m w + (j^ - yw , = m x ^ + f/ A ^ , 

where, p XfJjU v\ = as it obviously follows from (pUj) and (JSHJ)- 

4. Conservation laws and equations of motion 

The local conservation laws (fTTH) will then yield some restrictions on the coefficients p^ v 
and p XtiV [THHn]. First of all, the symmetry of t^ v implies that 

pV v = p u ^ p x ^ u = p Xu v 

Now, it is helpful to separate these coefficients in their components respectively parallel 
and orthogonal to the velocity v M : 

p^ = Mv»v u + p>*v v + p v v» + p^ , 

p x ^ = Q x v ^v u + Q^v" + Q Xv v» + Q Xflv , 
where all tensors and vectors other than are orthogonal to the velocity. The local 
conservation law (fT0|) then implies that 

<W + fyC + ^Ct + = • (42) 

Now, since 9^ t is the energy-momentum tensor of a free electromagnetic field, d^9^ t = 0. 
Similarly, the cross term contribution is 

<9,Cx = -FZJ» = ~ej drFZ(z) v,{r) S(x - z{r)) (43) 



(41) 
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and (see Appendix A [equation lJ77jl ] for details) 

2 



3^6% = -e 2 / dr [aV - d u ] 8(x - z{r)) . 



(44) 



Finally, using (jH| and after several integrations by parts, we also obtain 



dr 



J- (Mv u +p u + a A [Q V + Q Xu ]) 5(x - z) 



+ (v»p» + P f + r x ^[Q x v u + Q Xv ]j - z) 

+ (Q XfW + Q x »v») d XfM 5(x - z)] 
and, substituting (JIBJ), (jUJ) and (JH^) into l|35j> . we arrive at 
0= Jdr |J- (M/ + / + a A [QV + Q A "]) 



(45) 



+ -e 2 (aV-d") -F^j^-z) 
+ (Q A ^ + <9 v 5(x - z)] , 



(46) 



where F v = eF^ k (z)v li . 

As the derivatives of (^-functions in the r.h.s. are contracted with tensors that are 
transversal to the worldline, each term must vanish separately and therefore 

F v , (47) 



^- (Mv v + p v + a x [Q x v v + Q Xu }) + \e 2 {a 2 v u - a h 

(XT O 



d 



v^ + p^ + fj^[Q x v u + Q Xu ) 



= 0, (48) 
q(\h)v + q(M) v » = o . (49) 

Since Q x ^ and Q Xflv are orthogonal to v x and Q A/i!y = Q A ^, equation ijlHJ) implies 



that 

Q (v) = and Q A ^ = . 
Substituting this into iflHfl . we obtain 

= - Q M + v^Q x cb X - Q» x a x , 

= _ q>v _ qp> + v "Qi* ax + ^Q A -a A . 

Since p^ is symmetric and Q tJ,u is skewsymmetric, it follows that 



(50) 

(51) 
(52) 

(53) 
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and 

qh> = -q\m q ^ _ 2v [ ^Q u]x a x . (54) 

Finally, substituting (|5T)j) . (|5T|) and ijHIJ) into l|ITj) . after a short manipulation we arrive 
at 

^-([M+ 2Q x a x ] v v - Q u + 2Q Xu a x ^ + jje 2 (aV - a") = F u . (55) 

On the basis of solely the conservations of energy-momentum and angular 
momentum we have thus found that 

(a) the quantities M, Q x . . . , Q XfiU in equations (|4*T| can be written in terms of only 
ten independent particle variables: M, Q x and Q 3 that, 

(b) together with the worldline variables z m (t), f M (r), . . . are subject to the differential 
system (l5l)-(l55l). 

4-1. Total momentum and angular momentum 

Next, to have a clue of the physical meaning of M, Q x and Q Xu , we examine the total 
linear and angular momenta. 

The total linear momentum contained in the hypersurface T = {r = constant} in 
the optical coordinates (JBJ), i. e. the future light cone with vertex in z M (r), is 

P"(t) = I dEut^ with dZ u = -k„d 3 X. (56) 



jt 

Including now (|39|) . we have that the total momentum P M results from three 
contributions: 

F" = P» + P£ ix + P^ t , 

where P^ ix and P ext respectively come from the cross term 9^ ix and the external field 
term 9^ t in the energy-momentum tensor, and 

Pg = -Jd 3 XK{C + 0%) (57) 

is the contribution from the charge, i. e. the charge and its inseparable self-field. 

On substituting IjSSjl . ljS5|> . I^Tljl and JUJ into (|H7j> . after a little calculation we 
obtain 

P£ = Mr;' 1 - + ^(QV + g A ^) . (58) 
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Similarly, the total angular momentum in the hypersurface T, 

J^( T ) = - J d 3 Xk a {xH va - x"tr) , 

comes from three contributions as well: J^ v = J£ v + J^ x + J^ t . A similar calculation 
yields the point charge contribution 

TV — ^T> v — ? v -I- 

where 

= -2QK" ] - 2Q l ^ ] (59) 

is the particle internal angular momentum. The second term on the r. h. s. is 
orthogonal to the velocity and is the spin of the particle. On its turn, the possibility 
that Q M = —v v S^ u 7^ is related with the fact that the center of motion [TH] does not 
necessarily lies on the particle's worldline. 

To model a spinless charge, we choose Q^ v = 0. Equation (IH3I) then yields 

Q" = Qa" (60) 

and (JHKJ) can be further simplified to: 

d /r.. n „ • \ 2 



dr 



( [M + 2Q x a x ] v v - Q") + -e 2 (aV - a") = . (61) 

This agrees with the equation obtained by Honig and Szamosi: (ffiTf is equation (7) 
in [TT], with m = M + 2Qa 2 — Q, R = 2Q and S = Q. Lorentz-Dirac equation is a 
particular case for Q = 0. 

4-2. Summary 

A classical spinless point charge is therefore described by 
(a) the electric current density (l20l 

f = e [ drv^(T)S(x- z(t)), 



where the electric charge e is a constant scalar, and 
(b) the total energy-momentum tensor (|39"j) 

ir = c + K + Ct + Cx 
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where 9^ and t% u are respectively given by (j3l$| and (jlOJ, with 

jT = (M + 2Q A a A )*;'V - 2^- (Q^^) + Q'V . (63) 

ar 

The scalar variables M and g, together with the worldline z^ij) are subject to equation 
(jfiTJ) . which has been derived on the only basis that linear and angular momenta are 
conserved, supplemented with the point limit and the assumption that the particle is 
spinless. 

5. The equation of motion 

Equation (|6T|) does not yield the law of motion yet. Indeed, it consists of four equations 
for five unknowns, namely, M, Q and z M with the constraint v^v^ = — 1. The motion of 
the particle is therefore underdetermined. 

This should not be surprising. The problem in dynamics of continuous media for 
e > 0, as we have posed it, is itself underdetermined, because no constitutive equation 
has been assumed for the material sustaining the electric charge, contrary, for instance, 
to what is done in [HUES], were it is assumed that the charge is rigidly distributed over 
a spherical shell of radius e. 

Instead of advancing a matter constitutive equation for e > 0, then reexamining 
the problem and taking the limit e — > to determine a final equation of motion, we 
shall directly posit a constitutive relation connecting M, Q and the worldline invariants 
(curvature, torsion, etc.). 

Notice that, although it is the simplest choice and looks suitable for an elementary 
charge, a prescription like Q = is not an appropriate constitutive relation. Indeed, 
with a choice like this, (|6T|) becomes Lorentz-Dirac equation which leads to the dilemma 
of solutions that are either preaccelerated or runaway. 

We shall base our guess of a constitutive relation on the requirements that 

(a) it connects M, Q, a u and maybe some of their derivatives, 

(b) when a u , Q and also their derivatives vanish, then M = mo, and 

(c) if the point charge is acted by an external force F v that vanishes for r < and for 
t > Ti, then: 

• a v {r) = 0, M(r) = m and Q(t) = for r < and 
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• a u — > 0, M — > m and Qa u — > asymptotically in the future. 

(The proper mass has the same value mo in the infinite past and future, because we are 
assuming that the particle "identity" is finally preserved.) 

5.1. Rectilinear motion 

To see whether a constitutive relation can be prescribed so that (|6TJ) admits solutions 
that are neither runaway nor preaccelerated, we shall examine the case of rectilinear 
motion. (Recall that even in this simple case Lorentz-Dirac equation is not satisfactory.) 

Consider a point charge that initially is unaccelerated and free. Then, during the 
interval < r < T\, it is acted by an external force in a constant direction along the X 1 
axis. The charge worldline will remain in the plane X 1 X A in spacetime and therefore, 

dv^ , da* 1 . o ,. 

— — = a cr and — — = aar + a , 
dr dr 

where at 1 is the unit vector parallel to a M , i. e. the first normal to the worldline. The 
coefficients p^ v and p Xfll/ in equations (|62|) and (|63| . i. e. the particle's contribution to 
the energy-momentum tensor are 

p x ^ = qa x v ii v v , Q x = Qa x , (64) 

=M//-g(aV + aV)+gaa' 1 a F , (65) 

with q = Qa. 

In this case, the only non-vanishing components of equation (|6T|) are 



(II «") ^-(M + qa) = aq, 

2 

(l/) a (M + qa) — q e 2 d = F . 

3 



(66) 



These two equations must be supplemented with a constitutive relation M = 
M(a, q, q) in order that evolution is determined. The phase space is therefore 
coordinated by (a, q, q). 

We would expect that while the charge is not acted by any force, F(t) = 0, 
— oo < r < 0, then it remains in a state of uniform rectilinear motion and the energy- 
momentum tensor is the one corresponding to a free particle together with its Coulomb 
field, i. e. equations (|39l) . (j62|) and (|63|) with 

a(r) = 0, M(r)=m , q(r) = q{r) = , -oo < r < (67) 
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If an external force is then switched on: F(r) ^ 0, < r < r l5 then a, M, q and q 
evolve according to (|66|) with the initial data inferred from (|67]l and the continuity of 
the orbit in phase space. This determines 

a(r) , M(r) , q(r) and q(r) for < r < n (68) 

After that the particle is not acted by a force any more and what we would expect is 
that it asymptotically tends towards a free state, i. e. 

a(r) , M(t) m , g(r) -> , g(r) ^ for r^oo 

(with the same asymptotical value m for the mass, in order that the particle's "identity" 
is preserved). 

A way to achieve this behaviour consists in that the dynamical system (f66|) 
supplemented with the constitutive relation has only one equilibrium point for a = 
q = q = 0, which is asymptotically stable and M(0, 0, 0) = m . 

5.2. A dynamical system 

2e 2 

Using the constant r = , we introduce the new dimensionless variables 

3m 

r M + qa q 

t = — , 1 + fi = , a = ar , p = (69) 

t m m T 

and reduce (j66|) with F = to the simpler equivalent system 

jj! — ap' , p" + a' = a(l + p) , p, = p,(a, p, p') , 

where 'prime' means «derivative with respect to t». 

Then, by differentiating the constitutive relation and introducing the variable 
x = p' + a, we obtain 

p' = x — a , I 

x> = a(l + pi) , > (70) 
a' = A(a, p, x) , J 

where 

A(a,p,x) = - [(x- a) (a - pL p ) - ap, x (l + p)]. 

This dynamical system is already in normal form and is defined in the entire phase 
space provided that the function A(a, p, x) has no singularities. Particularly, if we choose 
so that is a solution of 

A (a, p, x) p, a + (x — a) p, p + a (1 + p,) p, x = a(x — a) (71) 
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with A (a,p,x) = la + pp + rx (Z, p and r constant) and p(0, 0, 0) = 0, then the 
dynamical system (|To"j) becomes 



cZ_ 
eft 



P 

\ x / 



^ I p r ' 

-1 1 
\ 1 0/ 



/a\ / \ 

p + 
\ x J \ pa J 



(72) 



If p ^ 0, the equilibrium points are 

Pj : a = p = x = , 
Pj/ : x = a = a , Po 



Z + r 
P 



«o and /i(ao, po, «o) = — 1 • 



Moreover, the constants I, p and r can be chosen so that the characteristic equation at 
Pi, 

X 3 - IX 2 + (p - r)X - p = , 

has three negative solutions and hence Pj is an asymptotically stable equilibrium point. 

In Appendix B [equation (|8TJ) ] we see how a solution p = p(a, p, x) of equation (fTTf 
that vanishes at Pj = (0, 0, 0) can be perturbatively obtained and is valid at least in a 
neigbourhood of this phase point. 

Now, (|69|) can be used to obtain the constitutive equation 



M = m — qam p I ar 



q q 

,ar H 

m T m 



(73) 



This, together with equations (fBfijl . determines a motion of the charge that is free of 
both preacceleration and runaways, provided that the force F acts only during a finite 
interval of time. Indeed, if the charge is unaccelerated in past infinity it remains so 
until its state is altered because F has started to act. Then, when the force ceases, the 
charge tends to the asymptotically stable equilibrium point a = 0, q = q = 0, at least if 
the system was close enough when the force dissapeared. 



6. Conclusion 

By studying the energy-momentum balance of a classical point charge with the 
electromagnetic field, we have obtained that 

(a) the total energy-momentum tensor consists of (i) a regular part, which comes from 
the external field contribution plus the regularization of the self-field contribution, 
and (ii) a singular part, with support on the charge worldline. 
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(b) This singular part depends on two scalar coefficients M(r) and Q(r) and on the 
worldline variables a^(r), . . . 

(c) These variables are constrained to fulfill the Honig-Szamosi equation [TT], i. e. (JHTJ). 

Lorentz-Dirac equation is obtained only if the constitutive relation Q = is set by 
hand. The well known troubles that suffers the Lorentz-Dirac equation are due to this 
bad choice rather than to energy-momentum conservation itself. 

We have then seen that, at least in the case of rectilinear motion, it is possible to 
find a constitutive relation M = M(a, Q, Q) which, together with equation (|6TJ) yields 
an equation of motion for the point charge that is free from both preacceleration and 
runaways. That is, if a charge is initially at rest, with proper mass mo, and is acted by 
an external force which lasts only a finite interval of time, then there is no acceleration 
before the force starts and, when its action ceases, the motion tends asymptotically to 
be rectilinear uniform and the proper mass tends to m . 
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Appendix A: Detailed computation of Eq. ()44R 

Using the definition (|33j) . we have that Wip G D'(IR 4 ) 



The first term vanishes because there is no current in p > e and, applying Gauss theorem, 
the second one yields 




(74) 



Since 6^(x) is summable for p > e, the first integral on the r.h.s. becomes 





(75) 
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where (an) = a x n x and d 2 Q is the solid angle element. Using then equation (J2nj) and the 
Taylor expansion |3] (p(z + ek) = <p(z) + ek x d x p(z) + ^e 2 k fJ, k x d fl \p(z) + 0(e 3 ), equation 
(1751) yields 

r°° r /72n f P 2 . 
J 1= J dr J — {-—(v"e(an)[p> + ek x d x p>] 

+ n u [l + t(an)} [<p + tk x d\p + ^e 2 k»k% x p]) 

+ j[a v ~ (a 2 )n u ] [ip + ek x d x <p] + e 2 [a 2 - (an) 2 }k u y> J + 0(e) . 
On integration with respect to d 2 fl and using that 

d 2 ttn u = [ d 2 {ln u n»n x = and [ d 2 ttn u n» = ^f)^ , 



we arrive at 



2 roc / 1 \ Op 2 f 00 

J! = - J dr (a«V - -rTd^J + — J dr [aV - a v ]p . (76) 

It is straightforward to check that the first term on the r.h.s. exactly compensates the 
second term on the r.h.s. in (17411 . Therefore we have 



dj)% = jje 2 J dr [aV - a v ] 8{x - z(r)) . (77) 

Appendix B: The constitutive relation 

We have to solve equation ijTTjl 

(la + pp + rx) p a + (x — a) p p + a (1 + p) p x = a(x — a) (78) 

with the "initial condition" p(0, 0, 0) = 0. 

It is easily seen that this equation admits a perturbative solution like 

oo 

71=1 

fj,( n ' being a polynomial in the variables a, q, x which is homogeneous and has degree 2n. 
If we write 

D = (la + pp + rx)d a + (x — a)d p + ad x 
then equation (T7K1) yields the hierarchy: 

Dp (1) = a(x - a) , (79) 

oo 

n > 1 Dp {n) = - J2 V {n ~ 8) a d ^ {s) ■ ( 80 ) 

s=l 
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The lowest order is relatively easy to solve and yields: 

fi = — — — - [(p — r)a 2 + p 2 p 2 + (r 2 + p + ri)x 2 — 2pax + 2rppx] +0(4). (81) 
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